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, ? . GRPC Semi-Digital HCAL is a solid option for the PFA oriented calorimetry of the 

International Linear Collider. Together with the hardware, the software developments 

On \ is progressing steadily. The stauts and plans for the GRPC SDHCAL software devel- 

opment are presented, as well the first order digitization module for the GRPC and the 
display program DRUID (Display Root module Used for ILD) have been introduced. 

^ . 1 Introduction 

f5 . The Semi-Digital Glass RFC technology is a solid option for the ILD detector, because 

the RFC is a well-understood technology with many advantages, such as high efficiency, 

homogeneous, robust and low cost; as a gaseous detector, the RFC is almost free of neutron 

hits; and most importantly, the semi-digital technology codes every channel with only two 

^^^. bits, allows a granularity as high as 1 channel per square cm with relatively low electronic 

^ r^ ■ cost. Good performance is expected in jet energy resolution. 

_ Ph! The GRFC SDHCAL (GRFC Semi-Digital HCAL) group is an international collabration 

consists of more than ten institutes. The aim is to have a good understanding to the GRFC 

SDHCAL detector and to proof the GRFC SDHCAL is a feasible option to the ILC detector. 

^ ■ In order to understand the detector performance and to optimize the experimental set- 

f""^ . ting, prototypes (see Fig. [1]) are constructed and tested with various options. We have the 

^i: I so called mini-DHCAL with Ik channel, the square meter with 10k channel, and the ANR 

^sf^ . cubic meter with 400k channel is under constrtuction. Lots of new technologies will be 

tested on the cubic meter, such as power pulsing, embedded electronics and self-supporting 

mechanics. Besides, the analysis of the test beam data will be used to vahdate the Geant 4 

^^ ■ hadronic simulation. 

Now let's start from an introudction on the status and plan of GRFC Semi-DHCAL 
softwares. 
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^ . 1.1 Staus and plan for the GRPC SDHCAL softwares 
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As well as for the hardware, the GRFC SDHCAL group has made great efforts on the 
software development. 

For the test beam experiment, a software chain (DAQ, reconstruction, analysis) based 
ROOT has been developed. The data format will be updated to LCIO standard, and 
currently people are focusing on the R & D of the new DAQ system. 

A Geant 4 simulation package with the cubic meter prototype geometry has been de- 
veloped, and a valid concept for the ILD with a la Videau DHCAL Geometry and GRFC 
sensor layer is currently availalbe in Mokka, the standard ILC detector simulation software. 



*On the behavior of SDHCAL collabration 
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Figure 1: GRPC SDHCAL Prototypes: mini-DHCAL, square meter and ANR cubic meter 

The digitization is important for the GRPC SDHCAL, as we beheved the semi-digital 
option could kept more useful information than the pure digital HCAL. A first order digiti- 
zation module with the cosmic ray experimental input has been developed, which converts 
the shower particle energy deposition in each cell to the quantity of induced charge on the 
electrinic nodes. Preiliminary result is optimistic. The digitization module has been in- 
tegrated into standard Marlin framework, and it will be upgraded with multiplicity and 
saturation correction. 

Because GRPC SDHCAL utilises totally different technology and detector geometry 
comparing to the AHCAL concept, its reconstruction algorithm need to be developed and 
optimized, which is regarded as one of the central task for SDHCAL software development. 
People are studying the hadronic shower development, aiming at develop an optimized re- 
construction algorithm for the GRPC SDHCAL. The idea is to integrate the SDHCAL 
reconstruction algorithm into the full Pandora reconstruction chain [l] , to achieve the best 
performance in jet energy resolution. 

A 3d display program DRUID based on root ^ TEve class ( developed for the LHC 
event display ) has been developed. It can visualize the information stored in standard ilc 
data file (slcio file) with various options and different detector geometries. With DRUID, 
better understanding of the ILD/test beam event and reconstruction algorithm performance 
could be achieved. Some examples are given in the following section. 

Central MC generation with particle gun events and ILC benchmark physics processes 
is also proposed, which will server as a basis for the algorithm development, performance 
analysis and geometry optimization study. 

Now let's fouse on the software tools of digitization and display. 

2 Digitization module for GRPC SDHCAL 

The digitization module is important for the GRPC SDHCAL option since people are eager 
to know if the SDHCAL option could successfully record more useful information than 
DHCAL. 

As shown in Fig. [21 a cosmic ray experiment has been performanced. Using the scin- 
tillator as trigger, the anaylog readout of one channel has been recorded, and reconstruced 
into the induced charge spectrum with about 2700 cosmic ray events. The induced charge 
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Figure 2: Cosmic ray experiment and Induced charge spectrum 



spectrum could be fit to a Polya function as signal and a Gaussian pedestal noise. The fitted 
parameters of the spectrum are used as the input for the digitization. 

The method of digitization is straightforward. First, readout the Monte Carlo truth 
information of energy deposition from simulation; second, estimate the number of ionizations 
that will be generated according to the energy deposition. The induced charge to one 
individual ionization is estimated from a Polya function (with the parameters calculated 
from the experimental input). And the total induced charge is assumed to be sum of the 
charge of all the ionization. 
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Figure 3: Correlation between induced charge and deopsited energy for lOOGeV Pion 

The correlation of induced charge and deposited energy is shown in Fig. [3l Because 
of large uncertaintly of the RPC avanache development, there is large smearing in charge 
inducing at low energy deposition region, but statistically nice correlation has been observed 
at high energy deposition region, where the induced charge could successfully represnent the 
energy deposition information. In other word, SDHCAL do keeps more energy information 
comparing to the pure DHCAL. 
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Figure 4: Left: Simulated ZH event at ILD detector Right: Reconstructed C ALICE test 
beam event 



The current digitization module has not taken into account the multiplicity and satura- 
tion effects. Those corrections will be included in a upgraded version. 

3 DRUID: Displaying Root module used for ILD 

Druid [31 is a 3d display program used to visualize information stored in standard ilc data 
file (slcio format). Currently it is specified to support the ILD geometries with gear file as 
geometry input file, and has been tested on simulated ILD event and CALICE test beam 
event (Fig. H]). According to different collections stored in the Icio file, Druid could display 
different objects in different combinations and styles. 

Druid needs a pre installed ROOT v5. 22.00, GEAR vO.12.0 and LCIO vl.ll (or higher). 
It can be accessed at IN2P3 svn server or my personal webpage 



http : //llr . in2p3 . f r/~ruan/ILDDisplay with a brief manual. 



3.1 Objects to be displayed 

3.1.1 Detector geometry 

Currently Druid support the ILD full detector geometry with a la Videau HCAL or TESLA 
HCAL, and the frame work for CALICE test beam geometry. Different subdetectors could 
be easily mounted and dismounted by interactivly clicking on the GUI. 

3.1.2 Event information 

Event information is organized into different groups: 

MCParticle: MCParticles are displayed as tracks or arrows. 

Calorimeter Hit: simulated/digitized calorimeter hits, displayed as cuboids with different 
size, orientations and color. 

Tracker Hit: simulated/digitized tracker hits, displayed as points. 

ReconstructedParticle: reconstructed particles are displayed as tracks with the calorime- 
ter hits assigned to this particle. 
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Figure 5: 40 GeV pion shower in the ILD Calorimetry 

Here MCParticle and reconstructed particle cohection are divided into subgroups accord- 
ing to their PID/energy, and the hits groups are divided into subgroups according to the 
subdetectors. 

Different groups of objects could be hidden/displayed. Druid will remember the status 
of display /hide for different objects from last event, and keep the status when navigated to 
a new event. By default, Druid will always display the new object. 

3.2 Display options 

Druid display could be easily zoomed, rotated and shifted by the mouse. With a self defined 
GUI pannel, lots of interactive actions are supported in DRUID. 

3.2.1 Individual objects 

TEve allows us to attach text information on individual objects, which could be readout by 
picking up the object with the mouse. 

The rotation center could be selected, together with the zoom option, you can have a 
very closed and clear view of event details. Fig. [5] shows the shower created by 40GeV 
Pion, with the Monte Carlo truth information of the secondary particles generated inside 
the hardonic shower. Where you can see a 2.04GeV Pion is generated in the interaction, 
and pass through the HCAL as a mip. 

3.2.2 Option for calorimeter hits 

Druid allows various of options for the size and color of calorimeter hits, which could be 
specified in the option pannel. For example, simulated calorimeter hit could be coloured 
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Figure 6: Tau jet r ^ i^ + tt^ + 7r+ with different color options: from left to right, energy, 
track ID and random index on track 
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Figure 7: Tau jet r ^ i^ + tt^ + 7r+ at MCTruth and Reconstruction level 



according to the hit energy, the PDG of MCParticle that passing through or at the vertex (as 
origin), or a randomly generated color. As for the reconstructed particle assigned calorimeter 
hit, it could be coloured to the PDG information of reconstructed particle. 

Fig. [6] shows how a r jet looks like with different color option (at MC truth level). 

3.3 Analysis reconstruction algorithm performance 

By comparing the Monte Carlo truth objects and the digitized/reconstructed objects, DRUID 
provides a platform to analysis the reconstruction software performance. 

Fig. [71 shows an example with the same r jet at MC truth and reconstruction level. The 
left one is the reconstructed particle tracks (using PandoraPFA) with assigned calorimeter 
hits, where you can see Pandora found two gamma and one pion; The middle one overlay 
the MCParticle information, thus we know that Pandora is right about the gammas and 
pions, also some secondary particles generated inside the shower and the neutrino could be 
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found. The right one overlay the simulated calorimeter hits: some simulated calorimeter 
hits are dropped during the digitization/reconstruction. You can easily tell if this dropping 
is expected or not by reading the text information attached to each of those hits. 

3.4 Conclusion 

As a solid option for the ILC, the software development of GRPC SDHCAL technology is 
progressing at steady pace with the hardware developments. 

The first order of digitization module has been developed, and the prehminary result 
shows that SDHCAL can record more information than the DHCAL, ensures better perfor- 
mance. 

The display package DRUID has been developed, with which people could get better 
understanding to the ILD/test beam event as well as to the performance of reconstruction 
algorithms. 

The digitization module will be upgraded with realistic physics efforts (saturation, effi- 
ciency and multiplicity) and the DRUID will be updated with new geometries. As a cen- 
tral task for Semi DHCAL software development, people are currently foucs on the design 
of reconstruction algorithm, which will be integrated into the PandoraPFA reconstruction 
framework. 

4 Bibliography 
References 

[1] |http : //www . hep . phy . cam . ac . uk/~ thomson/pandoraPFA | 

[2] |http : //root . cern . ch| 

[3] |http : //llr . in2p3 . f r/~ ruan/ILDDisplay/DruidNote . pdf | 



LCWS/ILC2010 



-Event Navigation 




l-SimealoHit Colour According to: -i 



PDG of Tracl^ 



1.0 vj 
Size oftlie Simuliits: 
Ratio to sidesize = 2cm 



pPFOCaloHit Colour According to: 

~3 



PDG of Track 



1.0 

Size of the PRO clustered hits: 
Ratio to sidesize = 2cm 



^Digitized Calorimeter Hit Size - 
P Scale the calorimeter hits 
according to the logarithm 
of hit energy 



I 1° 

Size of the hits: 

Ratio to sidesize = £cm 



^OHCAL Colormap 
Thresholds: 
1 0, 1 00 a 1 KeV: 



pGlobal SF for Digitized Hit 
Cell Color Scale; 



10 




^t 

Hid-e 



Eve I Files Options 
-Event Navigation 



<ik> 



^4t> 



Go to 
Evt 



pSimCaloHit Colour According to: -i 



PDG of Tracit 

Size of the Simuhits: 
Ratio to sidesize = 2cm 




P 



j-PFOCaloHit Colo ur Acc ording to: 
PDG of Trac k ^ 



Size of the PRO clustered hits: 
Ratio to sidesize = 2cm 



P 



t 



Digitized Calorimeter Hit Size- 
n Scale the calorimeter hits 

according to the logarithm 

of hit energy 

Size of the hits: 

Ratio to sidesize = 2cm 



DHCAL Colormap 
Thresholds: 
10. 100 & IkeV: 



^ 



Global SF for Digitized Hit 
Cell Color Scale: 







H 



Browser Eve: 



Eve 



Files Options 



[7: Wlndo"^ Manager 
"^ Vle-^ers 
7 Scenes 
S- -^ [5^ Geometry scene 
\^' _J [^ Trackers D 
S"^ i^ Calorimeters D 

i-£]p^ECain 
I [^ Origin marker ■ 
B-^ I? Event scene 
? Event 

■Qr MCPARTICLE ■ 
F"^l~ Calorimeter Hit D 
R-^^jy ReconstructedParticie D 
Fl- ^ R FandoraPFOs ■ 
[^ LowE ■ 
j7 Gamma D 
_J [7 Neutral Had D 
P.I ly Muons ■ 
_J[^ Pions ■ 
J [^ Electrons ■ 

I [^' Charged Kaons I 

C] [7 Protons ■ 
_J [7 Klong D 
J [7 Neutrons D 

I [7 Ion D 

■ Pin Sim Calorimeter Hit D 
■Q|~SimTrackerHit D 
■OP Tracker Hit D 



ilii 



Viewer 1 
j^Hide 




.Styfe 

P HCal UEveCompound] 



TEveElement 



Show: r Sejf F Children I I 




Command 



